Graphite and graphene as perfect spin filters 
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Based upon the observations (i) that their in-plane lattice constants match almost perfectly and 
(ii) that their electronic structures overlap in reciprocal space for one spin direction only, we predict 
perfect spin filtering for interfaces between graphite and (111) fee or (0001) hep Ni or Co. The 
spin filtering is quite insensitive to roughness and disorder. The formation of a chemical bond 
between graphite and the open d-shell transition metals that might complicate or even prevent spin 
injection into a single graphene sheet can be simply prevented by dusting Ni or Co with one or a 
few monolayers of Cu while still preserving the ideal spin injection property. 

PACS numbers: 72.25.-b,73.43.Qt,75.47.-m,81.05.Uw,85.75.-d 



The observation [TJ [2] of giant magnetoresistance 
(GMR) in systems where the transmission through inter- 
faces between normal and ferromagnetic metals (FM) is 
spin-dependent has driven a major effort to study spin fil- 
tering effects in other systems. An ideal spin filter would 
allow all carriers with one spin through but none with 
the other spin. Interfaces with half-metallic ferromag- 
nets (HMFs) [3] should have this property but progress 
in exploiting it has been slow because of the difficulty of 
making stoichiometric HMFs with the theoretically pre- 
dicted bulk properties and then making devices main- 
taining these properties at interfaces [4 . 

If the nonmagnetic metal (NM) is replaced by an in- 
sulator (I) or semiconductor (SC), spin filtering still oc- 
curs giving rise to tunneling magnetoresistance (TMR) 
in FM|I|FM magnetic tunnel junctions (MTJs) and spin- 
injection at FM|SC interfaces. If the spin-polarization of 
the ferromagnet is not complete, then the conductivity 
mismatch between metals and semiconductors or insu- 
lators has been identified as a serious obstacle to effi- 
cient spin injection [5]. It can be overcome if there is a 
large spin-dependent interface resistance but this is very 
sensitive to the detailed atomic structure and chemical 
composition of the interface. Knowledge of the interface 
structure is a necessary preliminary to analyzing spin fil- 
tering theoretically and progress has been severely ham- 
pered by the difficulty of experimentally characterizing 
FM|I and FM|SC interfaces. 

The situation improved with the confirmation of large 
values of TMR in tunnel barriers based upon crystalline 
MgO [6l [7] which had been predicted by detailed elec- 
tronic structure calculations [8j E] . While the record val- 
ues of TMR - in excess of 500% at low temperatures [T0| - 
are undoubtedly correlated with the crystallinity of MgO, 



the nature of this relationship is not trivial [IT]. The 
sensitivity of TMR (and spin injection) to details of the 
interface structure [I2j [13] make it difficult to close the 
quantitative gap between theory and experiment. In view 
of the reactivity of the open-shell transition metal (TM) 
ferromagnet s Fe, Co and Ni with typical semiconductors 
and insulators, preparing interfaces where disorder does 
not dominate the spin filtering properties remains a chal- 
lenge. With this in mind, we wish to draw attention to a 
quite different material system which should be intrin- 
sically ordered, for which an unambiguous theoretical 
prediction of perfect spin filtering can be made in the 
absence of disorder, and which is much less sensitive to 
interface roughness and alloy disorder than TMR or spin 



TABLE I: Lattice constants of Co, Ni, Cu, and graphene. 
ahex = afcc/V 7 !. Equilibrium separation do for a layer of 
graphene on top of graphite, Co, Ni or Cu calculated within 
the local density approximation (LDA) of density functional 
theory with a = 2.46 A. The binding energy AE(d ) = E(d = 
oo) — E(do) is the energy (per interface unit cell) required to 
remove a single graphene layer from a graphite stack or from 
a Co, Ni or Cu (111) surface. W is the workfunction. 
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FIG. 1: Fee Fermi surface (FS) projections onto a plane per- 
pendicular to the [111] direction for Co majority (a) and mi- 
nority (b) spins, for Ni majority (c) and minority (d) spins 
and for Cu (e). The number of FS sheets is shown by the 
colour bar on the right. For graphene and graphite, surfaces 
of constant energy are centred on the K point (f). 



injection. 

We begin by observing that the in-plane lattice con- 
stants of graphene and graphite match the surface lattice 
constants of (111) Co, Ni and Cu almost perfectly. From 
Table [TJ it can be seen that Ni is particularly suitable 
with a lattice mismatch of only 1.3%. The second point 
to note is that the only electronic states at or close to 
the Fermi energy in graphene or graphite are to be found 
near to the high symmetry K point in reciprocal space 
where Co and Ni have states with minority spin charac- 
ter only. The absence of majority spin states in a large 
region about the K point is made clear in the (111) Fermi 
surface (FS) projections shown in Fig.jl] The (0001) FS 
projections for hep Co are qualitatively the same. It fol- 
lows that in the absence of symmetry-lowering (resulting 
from disorder, interface reconstruction etc.) perfect spin 
filtering should occur for graphite on top of a flat Ni or 
Co (111) surface. 

The effectiveness of the spin filtering is tested for a 
current-perpendicular-to-the-plane (CPP) structure with 
n graphene layers sandwiched between semi-infinite Ni 
electrodes. The spin dependent transmission through 
this Ni|Gr n |Ni junction is calculated using a first- 
principles tight-binding muffin tin orbital (TB-MTO) 
wave-function matching scheme [TT] [18] for parallel (P) 
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FIG. 2: Conductances Gf in (v), Gp aj (A), and G a AP (x) of 
a Ni|Gr n |Ni junction as a function of the number of graphene 
layers n for ideal junctions. Inset: magnetoresistance as 
a function of n for: (circles) ideal junctions; (diamonds) 
Ni|Gr n |Cu5oNi5o|Ni junctions where the surface layer is a dis- 
ordered alloy; (squares) Ni|Gr n |Ni junctions where the top 
layer of one of the electrodes is rough with only half of the 
top layer sites occupied (sketch). Results for roughness and 
disorder are modelled in 5 x 5 lateral supercells and averaged 
over 20 disorder configurations as described in Refs 1 131 1181 
For the rough surface layer, the error bars indicate the spread 
of MR obtained for different configurations. The supercell 
conductances are normalized to the lxl surface unit cell 
used for the ideal case. 



and antiparallel (AP) orientations of the Ni magneti- 
zations. The atomic sphere (AS) potentials are calcu- 
lated self-consistently within density functional theory 
for atomic structures determined by total energy mini- 
mization (see below). The conductances G P and G G AP 
are shown in Fig. [2] for the minority and majority spin 
channels, a = min, maj. Gp aj and G a AP are strongly 
attenuated while, apart from an even-odd oscillation, 
Gp m is independent of n. The magnetoresistance MR 
= (Rap - Rp)/Rap = (G P - G A p)/G P rapidly ap- 
proaches 100%; see inset. We use the pessimistic defi- 
nition of MR because Gap vanishes for large n; usually 
the optimistic version is quoted [6j [9] [10] . Similar 
results are obtained for Ni|Gr n |Co and Co|Gr n |Co junc- 
tions. 

We envisage a procedure in which thin graphite lay- 
ers are prepared by micromechanical cleavage of bulk 
graphite onto a Si02 covered Si wafer [19] into which 
TM (Ni or Co) electrodes have been embedded and lay- 
ers of graphene are peeled away until the desired value 
of n is reached. Assuming it will be possible to realize 
one essentially perfect interface in a CPP geometry, we 
studied the effect of roughness and disorder at the other 
interface on MR (inset Fig.[2|. Replacing the top Ni layer 
with a NisoCuso random alloy only reduces the MR to 
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FIG. 3: Total energy E of a graphene sheet on fee Co, Ni 
and Cu (111) surfaces as a function of the separation of the 
graphene sheet from the top layer of the metal. The lowest 
energy "AC" configuration is sketched on the bottom right. 

90% (900% in the optimistic definition). Extreme rough- 
ness, whereby half of the Ni interface layer is removed 
at random, only reduces the MR to 70%. The momen- 
tum transfer induced by the scattering is apparently in- 
sufficient to bridge the large gap about the K point in 
the majority spin FS projections. Alternatively, it may 
be possible to prepare two separate, near-perfect TM|Gr 
interfaces and join them using a method analogous to 
vacuum bonding [20 . 

Graphite has a large c-axis resistivity [21]. If one of the 
TM|Gr interfaces is ideal and the graphite layer is suffi- 
ciently thick, then it should be possible to achieve 100% 
spin accumulation in a high resistivity material making 
it suitable for injecting spins into semiconductors [5 . Be- 
cause carbon is so light, spin- flip scattering arising from 
spin-orbit interaction should be negligible. 

The results shown in Fig. [2] were calculated for the 
lowest energy "AC" configuration of graphene on Ni cor- 
responding to one carbon atom above a Ni atom (the 
surface "A" sites) while the other is above a third layer 
Ni "C" site. A and C refer to the conventional ABC 
stacking of the layers in an fee crystal (AB for an hep 
structure). The CPP spin filtering should not depend 
on the details of how graphite bonds to the metal sur- 
face as long as the translational symmetry parallel to the 
interfaces is preserved. This is confirmed by explicit cal- 
culation for the "AB" and "BC" bonding configurations 
for varying graphene-metal surface separation d. 

The nature of the bonding may well play an impor- 
tant role in realizing such an interface experimentally. 
In Fig. [3] we show the total energy of a graphene sheet 
on TM = Co, Ni and Cu (111) surfaces as a function 
of d where the zero of energy has been chosen so that 
E{d = oo) = for an uncharged graphene sheet. The 
density functional theory (DFT) calculations were car- 
ried out using the projector augmented wave (PAW) 
method [22j [23] , a plane wave basis set and the LDA, 
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FIG. 4: Majority and minority spin band structures (green) 
of a single graphene layer absorbed upon (both sides of) a 
13 layer (111) Ni slab for a BC configuration with d — 3.3 
A, and an AC configuration with d — 2.0 A. The bands re- 
plotted in black using the carbon p z character as a weighting 
factor are superimposed. The Fermi energy is indicated by 
the horizontal dashed line. 



as implemented in the VASP program [24j [25] . Super- 
cells containing a slab of at least six layers of metal atoms 
with a graphene sheet adsorbed on one side of the slab 
and a vacuum region of - 12 A were used. The Brillouin 
zone of the (lxl) surface unit cell was sampled using a 
36 x 36 k-point grid. The plane wave kinetic energy cut- 
off was 400 eV. To avoid interactions between periodic 
images of the slab a dipole correction was applied [26 . 
The atoms in the metal layers were fixed at their bulk 
positions. The experimental lattice constant of graphene 
a — 2.46 A is used as the lattice parameter ah ex for Co, Ni 
and Cu. To plot the band structures (Fig. [4]), a 13 layer 
slab with graphene absorbed on both sides was used. 

The most prominent feature of Fig. [3] is the prediction 
of a weak minimum in the binding energy curve for Cu 
of about 0.07 eV at an equilibrium separation d ~ 3.2 
A and deeper minima of 0.37 and 0.32 eV respectively 
for Co and Ni, at a smaller equilibrium separation of 
d ~ 2.0 A. In agreement with a recent first-principles cal- 
culation [27] and experiment [I5j [28] for graphene on Ni, 
we find that the lowest energy corresponds to an AC con- 
figuration. The finer details of the total energy surfaces 
depend on the choice of exchange-correlation potential, 
relaxation of the metal substrate, choice of in-plane lat- 
tice constant etc. and will be presented elsewhere. We 
restrict ourselves here to properties which do not depend 
on these details. 

The electronic structure of a single graphene layer does 
depend on d. For the less strongly bound BC configura- 
tion of Gr on Ni, the equilibrium separation is d ~ 3.3 
A and the characteristic band structure of an isolated 
graphene sheet is clearly recognizable; see Fig. [4] For the 



4 



lowest energy AC configuration, the interaction between 
the graphene sheet and Ni surface is much stronger, a gap 
is opened in the graphene derived p z bands and there are 
no graphene states at the K point in reciprocal space at 
the Fermi energy for the minority spin channel. This may 
prevent efficient spin injection into graphene in lateral, 
current-in-plane (CIP) devices [29]. However, there is a 
simple remedy. If a monolayer (or several layers) of Cu 
is deposited on Ni, graphene will form only a weak bond 
with Cu and the Fermi energy graphene states at the K 
point will only be weakly perturbed. Cu will attenuate 
the conductance of both spin channels because Cu has no 
states at the K point but will not change the spin injec- 
tion properties as long as it is sufficiently ordered as to 
preserve the translational symmetry; a completely mixed 
CuNi top layer reduces the MR in a Ni|Gr n |Cu5oNi5o|Ni 
junction only slightly (Fig. Cu will also oxidize less 
readily than the more reactive Ni or Co. Spin-flip scat- 
tering in a thin layer of Cu can be neglected. The weaker 
bonding of graphene to Cu may also have practical ad- 
vantages in sample preparation. 

Finally, we remark that graphene may exhibit curious 
bonding properties to a (111) surface of a Cui-^Ni^ or 
Cui-^Co^ alloy; as a function of increasing concentra- 
tion x, the weak minimum at d ~ 3.2 A will evolve into 
a deeper minimum at d ~ 2.0 A with the possibility of a 
double minimum occuring for some range of concentra- 
tion x\ a propensity to form a second minimum is already 
evident in the binding energy curve for Cu. Calculations 
are underway to examine this possibility. 

Planar interfaces between graphene and close-packed 
Co, Ni, or Cu represent a very flexible system for study- 
ing the influence of atomic and electronic structure on 
electrical contact with graphene related systems such as 
carbon nanotubes [30] where the nanotube geometry is 
very difficult to model using materials specific calcula- 
tions [31 . The binding energy curves in Fig. [3] and elec- 
tronic structures in Fig. [4] show that the closer proxim- 
ity resulting from stronger bonding does not necessarily 
lead to better electrical contact if bonding removes the 
carbon-related conducting states from the Fermi energy. 

Motivated by the recent progress in preparing and 
manipulating discrete, essentially atomically perfect 
graphene layers, we have used parameter- free, materi- 
als specific electronic structure calculations to explore 
the bonding and spin transport properties of a novel 
TM|Gr n system. We predict perfect spin filtering for 
ideal TM|Gr n |TM junctions with TM = Co or Ni. 
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